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Abstract. It is known that unstable periodic orbits of a given 
map give information about the natural measure of a chaotic at- 
tractor. In this work we show how these orbits can be used to 
CO ' calculate the density function of the first Poincare returns. The 

^ I close relation between periodic orbits and the Poincare returns al- 

lows for estimates of relevant quantities in dynamical systems, as 
the Kolmogorov- Sinai entropy, in terms of this density function, 
■^j^ ■ Since return times can be trivially observed and measured, our 

I approach to calculate this entropy is highly oriented to the treat- 

(<~^ i ment of experimental systems. We also develop a method for the 

^\ ' numerical computation of unstable periodic orbits. 

o ■ 



, 1. Introduction 

>< 

■ Knowing how often a dynamical system returns to some place in 

^ phase space is fundamental to understand dynamics. There is a well 

established way to quantify that: the first Poincare return (FPR), 
which measures how much time a trajectory of a dynamical system 
takes to make two consecutive returns to a given region. Due to their 
stochastic behaviour, given a return time it is not feasible to predict 
the future return times and for that reason one is usually interested in 
calculating the frequency with which the Poincare returns happen, the 
density of the first Poncare returns (DFP). 

This work explains the existence of a strong relationship between un- 
stable periodic orbits (UPOs) and the first Poincare returns in chaotic 
attractors. Unstable orbits and first Poincare returns have been usu- 
ally employed as a tool to analyse and characterise dynamical systems. 
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With our novel approach we can calculate how frequently returns hap- 
pen by knowing only a few unstable periodic orbits. Additionally, such 
relation allows us to easily estimate other fundamental quantities of 
dynamical systems such as the Kolmogorov-Sinai entropy. 

Our motivation to search for a theoretical and simple way of calcu- 
lating the distribution of Poincare return times comes from the fact 
that they can be simply and quickly accessible in experiments and also 
due to the wide range of complex systems that can be characterized 
by such a distribution. Among many examples, in Ref. [1] the return 
times were used to characterize a experimental chaotic laser, in Refs. 
[21 E] they were used to characterize extreme events, in Refs. [H [S] they 
were used to characterize fluctuations in fusion plasmas, and in Ref. 
[0] a series of application to complex data analysis were described. 

In addition, relevant quantifiers of low- dimensional chaotic systems 
may be obtained by the statistical properties of the FPR such as the 
dimensions and Lyapunov exponents [H [8] and the extreme value laws 
[9]. For most of the rigorous results concerning the FPR, in particular 
the form of the DFP p^, one needs to consider very long returns 
to arbitrarily small regions in phase space, a condition that imposes 
limitations into the real application to data sets. 

We first show how the DFP can be calculated from only a few UPOs 
inside a finite region. Then, we explain how the DFP can be used to 
calculate quantities as the Kolmogorov-Sinai entropy, even when only 
short return times are measured in finite regions of the phase space. 

Our work is organized as follows: We first introduce the work of 
Ref. [11], which relates the natural measure of a chaotic attractor to 
the UPOs embedded in a chaotic attractor. The measure of a chaotic 
attractor refers to the frequency of visits that a trajectory makes to a 
portion of the phase space. This measure is called natural when it is 
invariant for typical initial conditions. This appears in Sec. [21 along 
with the relevant definitions. In Sec. |3] we define /9(r, S) the density 
of first Poincare returns for a time r to a subset S of phase space and 
we study the relation between the UPOs and this function. This can 
be better understood if we classify the UPOs inside S as recurrent and 
non- recurrent. Recurrent are those UPOs that return more than once 
to the subset S before completing its cycle. Non-recurrent are UPOs 
that visits the subset S only once in a period. While in the calculation 
of the natural measure of S one should consider the two types of UPOs 
with a given large period inside it, for the calculation of the DFP for 
a time r one should consider only non-recurrent UPOs with a period 
r. Sec. m is mostly dedicated to show how to calculate p(r, S) even 
when not all non-recurrent UPOs of a large period are known. Such a 
situation typically arises when the time r is large. We have numerically 
shown that the error of our estimation becomes smaller, the longer the 
period of the UPOs and the larger the number of UPOs considered. 
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Throughout the paper we illustrate results by presenting the calcu- 
lations for the tent map. Finally, in Sec. E] we show numerical results 
on the logistic map that support our approach. In particular, we ob- 
tain numerical estimates of the Kolmogorov-Sinai entropy, the most 
successful invariant in dynamics, so far. The estimates are obtained 
considering the density of only short first return times, as discussed 
in Sec. [51 The UPOs of period p are computed numerically as stable 
periodic orbits of a system of p coupled cells, a method described in 

2. Definitions and results 

Consider a d-dimensional map of the form Xn+i = F{xn), where 
X & Q <Z and Q represents the phase space of the system. Consider 
A G Q to represent a chaotic attractor. By chaotic attractor we mean 
an attractor that has at least one positive Lyapunov exponent. 

For a subset 5* of the phase space and an initial condition xq in the 
basin of attraction of A, we define h{xq, S) as the fraction of time the 
trajectory originating at Xq spends in S in the limit that the length of 
the trajectory goes to infinity. So, 

(1) ,(,„.s,= ,taH£M£Amii!l. 

n— >-oo n 

Definition 2.1. // fi{xo, S) has the same value for almost every Xq 
(with respect to the Lebesgue measure) in the basin of attraction of A, 
then we call the value fi{S) the natural measure of S . 



For now we assume that our chaotic attractor A has always a natural 
measure associated to it, normalized to have fi{A) = 1. In particular 
this means that the attractor is ergodic[TT]. 

We also assume that the chaotic attractor A is mixing: given two 
subsets, Bi and B2, in A, we have: 

lim n{B^ n F-^iB^)) = /x(5i)/x(52). 

n— >oo 

In addition, we consider A to be a hyperbolic set. 

The eigenvalues of the Jacobian matrix of the n-th iterate, F^, at the 
jth fixed point xj of F" are denoted by Ai^, X2j, ^uj, \u+i)j, ^dj, 
where we order the eigenvalues from the biggest, in magnitude, to the 
lowest and the number of the unstable eigenvalues is u. Let Lj{n) be 
the product of absolute values of the unstable eigenvalues at Xj. 

Then it was proved by Bowen in 1972 [12] and also by Grebogi, Ott 
and Yorke in 1988 |TT] the following: 

Theorem 2.1. For mixing hyperbolic chaotic attractors, the natural 
probability measure of some closed subset S of the d-dimensional phase 
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space IS 



(2) f^{S)= 



, n] 



where the summation is taken over all the fixed points Xj G S of F"'. 

This formula is the representation of the natural measure in terms 
of the periodic orbits embedded in the chaotic attractor. To illustrate 
how it works let us take a simple example like the tent map: 

Example 2.1. Let us consider F : [0, 1] — ?■ [0, 1] such that 

_ r 2x, 2/ X G [0, 1/2] 
^■^^ \2-2x, if X G]l/2,1] 

For this map there is only one unstable direction. Since the absolute 
value of the derivative is constant in [0, 1] we have Ljij) = L{t) = 2"^. 

For the tent map, periodic points are uniformly distributed in [0, 1]. 
Using this fact together with some of the ideas of G.H. Gunaratne and 
I. Procaccia [I3], it is reasonable to write the natural measure of a 
subset S of [0, 1] as: 

where N{t, S) is the number of fixed points of F'^ in S and N{t) is the 
number of fixed points of F'^ in all space [0, 1]. For this particular case 
we have N{t) = L{t) = Lj{T) and so 

yu(S') = lim — ^ ' = lim — ^ / = hm — --r- 

r->oo N[T) T^oo L{t) t^oo -^f— ^ Lj[T) 

and we obtain the Grebogi, Ott and Yorke formula. 



3. Density of first returns and UPOs 

In this section we relate the DFP, p(r, S), and the UPOs of a chaotic 
attractor. We show in Eq. ffTOl) that p(r, S) can also be calculated 
in terms of the UPOs but one should consider in Eq. ([2]) only the 
non-recurrent ones. 



3.1. First Poincare returns. Consider a map F that generates a 
chaotic attractor A G fl, where fl is the phase space. The first Poincare 
return for a given subset S G Q such that SnA 7^ is defined as follows. 

Definition 3.1. A natural number t, t > 0, is the first Poincare 

return to S of a point xq & S if F'^{xo) G S and there is no other 
T* < T such that F^*(xo) G S. 



DENSITY OF FIRST RETURNS, PERIODIC ORBITS AND KS ENTROPY 5 



A trajectory generates an infinite sequence, Ti, r2, Xj, of first re- 
turns where ti = r and is the first Poincare return of F"''(xo) with 



Ei—l 
n=l 



The subset S' of points in 5* C that produce FPRs of length r to 
S is given by 

(4) s' = s'{t, s) = {F~^{S) ns) - \J {F-^{S) n s) . 

0<j<T 

3.2. Density function. In this work, we are concerned with systems 
for which the DFP decreases exponentially as the length of the return 
time goes to infinity. Such systems have mixing properties and as a 
consequence we expect to find p(r, 5*) ~ /i(S')(l — fi{S)y'~^, where 
(1 — iJ,{S)y~^ represents the probabihty of a trajectory remaining r — 1 
iterations out of the subset 5*. We are interested in systems for which 
the decay of p(r) is exponential, i.e., p(r) oc e"""^. 

The usual way of defining p(r, 5), for a given subset S C Q, is by 
measuring the fraction of returns to S that happen with a given length 
r with respect to all other possible first returns [see Eq. (127|) ]. It is 
usually required for a density that 



J p{t, S)dT = 1. 



In this work, we also adopt a more appropriate definition for p(r, S) 
in terms of the natural measure. We define the function p(r, S) as the 
natural measure of the set of orbits that makes a first return r to S* 
divided by the natural measure in 5*. More rigorously, we have: 

Definition 3.2. The density function of the first Poincare return r 
for a particular subset S C Q such that p(S') is defined as 

(5) ^(-«)-^. 

where S' = S'{t,S) C S is the subset of points that produce FPRs of 
length t defined in Eq. 

Even for a simple dynamical system as the tent map, the analytical 
calculation of p(r, S) is not trivial. However, an upper bound for this 
function can be easily derived as in the following example: 

Example 3.1. Consider the tent map defined in example \2.1\ for which 
the natural measure coincides with the Lebesgue measure X, and let 
S C [0, 1] be a non-trivial closed interval. 

To have a return to S we only need to know the natural number 
n* such that F^* (S) = [0,1]. Since F is an expansion, this natural 
number always exists. To find it when X{S) = e > 0, we first solve 
the equation 2^* = 1/e and get x* = — log(e)/log(2), so we take n* = 
[— log(e)/log(2)] + 1, where [x] represents the integer part of x. Then 
n* is an upper bound for T^m, the shortest first return to S. 
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Most intervals S of small measure have large values of and 
Tmin ^ n* is a good approximation. A sharper upper bound for Tmin 
in S is the lowest period of an UPO that visits it. 

The set D = F'^* {S) fl S* 7^ represents the fraction of points in S 
that return to S (not necessarily first return) after n* iterations. Using 
Eq. ^ and since S' <Z D we have 

It is natural to expect that for r of the order of n* and close to Tmin 
we have p{r, S) < 2^^ . 

We can write this equation as p(r, 5*) < e^~'^^°s^^''^ = e'-"'^'^^-', where 
Ai = log(2) is the Lyapunov exponent for the tent map. In fact, in 
1991, G. M. Zaslavsky and M. K. Tippett 

[Ilj[T5] presented one formula for the exact value of 

p{t, S) . That result can only be valid under the same conditions that 
we have used previously, i.e. t ^ Tmin and for most sets of sufficiently 
small measure e, so that Tmin ~ n* . 

3.3. Density function in terms of recurrent and non-recurrent 

UPOs. Since our chaotic attractor A is mixing, the natural measure 
associated with A satisfies, for any subset S of nonzero measure: 

^^^^ = }Too ]I(s) ■ 

We can write the right hand side of the last equation, for any positive 
r, in two terms: 

with S' as defined in Eq. (jl]) and where S* = S*{S,t) is the set of 
points in S that are mapped to S after r iterations but for which r is 
not the FPR to S, so S' U S* = {S n F-^{S)) and S' n S* = 0. 

An UPO of period r is recurrent with respect to a set S* C if there 
is a point Xq € S" in the UPO with F"'{xq) G S* for < n < r. In other 
words, its FPR is less than its period. Thus, the UPOs in the set S* 
are all recurrent. We refer to them as the recurrent UPOs inside S. 

Associated with the recurrent UPOs in 5* we define 

and associated with the non-recurrent UPOs in S we define 
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where L^ij) and L^^{t) refer, respectively, to the product of the ab- 
solute values of the unstable eigenvalues of recurrent and non-recurrent 
UPOs of period r that visit S. 
Notice that, if /i(S') 7^ 0, 

hm = lim Hr{t,S) 

and 

(9) lim = fJ'NR{r,S) 

r-i>oo fi(^^ j r-i>oo 

since fj,{S*) / fj,{S) measures the frequency with which chaotic trajecto- 
ries that are associated with the recurrent UPOs visit 5* and fi{S')/ fi{S) 
measures the frequency with which chaotic trajectories that are asso- 
ciated with the non-recurrent UPOs visit 5*. 

Comparing Eqs. ([5]), and ([9]) we obtain the following: 

Main Idea: For a chaotic attractor A generated by a mixing uni- 
formly hyperbolic map F , for a small subset S <Z A, generated by a 
Markov partition and such that the measure in S is provided by the 
UPOs inside it, we have that 

(10) p(r, 5) ^ /iArij(r, 5), 
for a sufficiently large t. Moreover, 

IJ,{S) = lim [p(r, S) + ij,r{t, S)]. 



A Markov partition is a very special splitting of the phase space. 
For the purpose of better justifying Eq. (fTOj) . if a region C(r) belongs 
to a Markov partition of order r, then there is a sub-interval C(r) of 
C{t) that after r iterations is mapped exactly over C(r). Moreover, 
points inside C{t) make first returns to C{t) after r iterations. Then, 
/iij,(r, C(r))=0. As a consequence, for sufficiently large r we can write 
that fi[C{T)] p[t,C{t)]. 

But approximation (11 01) remains valid for a small nonzero r. The 
reason for that is the following: Notice that from the way Kac's lemma 
is derived (see Sec. 18. ip . Eq. ([2]) can be written as 

r°° p(T,S)dT 

< r > 

where < r > represents the average of the FPRs inside 5", since 
p(r, S)dT = 1. This equation illustrates that any possible existing 
error in the calculation of fi{S) by Eq. ([2]) is a summation over all er- 
rors coming from p(r, S) for all values of r that we are considering. As 
shown in Ref. [H], fi{S) can be calculated by Eq. ([2]) using UPOs with 
a small and finite period p. This period is of the order of the time that 
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the Perron-Frobenius operator converges and thus hnearization around 
UPOs can be used to calculate the measure associated with them. As 
a consequence, if fi{S) can be well estimated for p ^ 30 then p(r, S) 
can be well estimated for r << p. As we will observe, considering r 
small, of the order of 5, we get a very good estimation for p(r, S). 

In addition, we observe in our numerical simulation that 5* does not 
need to be a cell in a Markov partition but just a small region located 
in an arbitrary location in Q. 

We say that an UPO has FPRs associated with it if the UPO is non- 
recurrent. See that for every UPO there is a neighborhood containing 
no other UPO with the same period. If the UPO is non-recurrent then 
all points inside a smaller neighborhood will produce FPRs associated 
with this UPO in the sense that their FPR coincides with the UPO's. 
Consider Tmin as the shortest first return in S. 

Case r < 2rmm 

UPOs of period r are non- recurrent. This is illustrated in Fig. [1] 
(A), where Tmin = 7, for the logistic map (c = 4). In that picture we 
observe that for r < 14 all FPRs are associated with UPOs. Because 
of this fact /i(>S'*) = and then all the chaotic trajectories that return 
to 5* are associated with non-recurrent UPOs. So, p(r, S) ~ ^{S) and 
thus, p(r, S) ^ fiNR{r, S). 

Case T > iTrain 

We can have recurrent UPOs of period r, that do not have first re- 
turns associated with them. As a consequence ^(5*) > and recurrent 
UPOs contribute to the measure of S. This is illustrated in Fig. [1] (B), 
when r = 16. 

4. How TO CALCULATE THE DENSITY OF FIRST POINCARE 

RETURNS 

A practical issue is how to calculate UNnij^ S). There are two rele- 
vant cases: All UPOs can be calculated; only a few can be calculated. 

Assuming r to be sufficiently small such that all UPOs of period r 
can be calculated and sufficiently large so that Eq. ( TTU]) is reasonably 
valid, fJ^NRiTyS) can be exactly calculated and we can easily estimate 
p(r, 5") from Eq.([Tn]), using p(r, S) ^ fJ,NR{r, S). 

When r is large then, typically, only a few UPOs can be calculated. 
For this case, it is difficult to use Eq. ffTOjl to estimate p(r, 5*) since there 
will be too many UPOs. In order to calculate p(r, S) using fiNRi^, S) 
we do the following. First notice that 



(11) 



fi{S) = lim (p7Vij(r, S) + iir{t, S)). 

T— >00 
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Figure 1. This picture shows some UPOs inside S C 
[0, 1] and first Poincare returns for the logistic map, 
[xn+i = 4x„,(l — Xn)]- In this example = 7. For 
r < 14 all UPOs have FPRs associated with them. For 
r > 14 (as in (B) for r = 16) some UPOs are recurrent. 
Picture (B) is a zoom of picture (A). 

Considering then r sufficiently large we have that 

fi{S) ^ /i7Vij(r, S) + fiR{T, S) 

which can be rewritten (using Eq. ffTOl) which says that p(r, S) 
fJ'NRi'T, S), for finite r] as 

i^niT, s y 

This equation allows us to reproduce, approximately, the function 
p{t,S), for any sufficiently large r, only using the estimated value of 
the quotient 

fi{S) 

that is easy to obtain numerically, since not all UPOs should be calcu- 
lated but just a few ones with period r. We discuss this in 14. II below. 

4.1. How can we estimate fiji{T, S)/ fi{S)7 Considering a subset S 
and fixing r, we calculate a number t of different UPOs with period 
r (say, t = 50) inside S (It is explained in Sec. 16.51 how to calculate 
numerically UPOs with any period of a given map). These UPOs are 



(12) pir,S)^fiiS)-fxniT,S) = fxiS){l 
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calculated from randomly selected symbolic sequences for which the 
generated UPOs visit S. See that, for example, in the tent map, for 
r = 10 and S = [0, |], we may have 2^^/ 8 UPOs inside S and so, here 
50 UPOs inside S is, in fact, a very small number of UPOs. 

Now, we separate all the t UPOs that visit 5* into recurrent and 
non-recurrent ones and suppose that we have r recurrent and nr non- 
recurrent such that r + nr = t. So, r and nr depend on t and 5*. With 
these particular r(t, S) recurrent UPOs we use Eq. ([7]) and we obtain 



where L^ij) represents the product of the absolute values of the un- 
stable eigenvalues of the jth recurrent UPO within the set of r{t,S) 
recurrent UPOs. See that this quantity is not equal to fiR^r, S) since we 
are not considering all recurrent UPOs inside S but just a small number 
r{t, S) of them. We do the same thing with the nr{t, S) non-recurrent 
UPOs and obtain the quantity fiNRlT, S,nr{t, S)]. 

Finally, we observe that, for a sufficiently large t, we have 



where /i(r, 5, t) = IIrIt, S,r{t, S)]+ (InrIt, S,nr{t, S)]. Therefore, with 
only a few UPOs inside S we calculate an estimated value for p(r, S). 
This estimation is represented by pm and is given by 



Notice that, for a large r we will have more recurrent UPOs than 
non-recurrent ones and therefore the larger r is, the larger is the con- 
tribution of the recurrent UPOs to the measure inside S. 

4.2. Error in the estimation. To study how much our estimation in 
Eq. f ll3p depends on the number t of UPOs, we first assume that if all 
UPOs are known, the calculated distribution in Eq. (ITU]) is "exact", or 
in other words it has a neglectable error as when compared to the real 
distribution provided by Eq. (|5]). 

Then, the error in Eq. f|T3|) will depend on the deviation of the 
quotient 




r{t,S) 



PrIt, S,r{t,S)] _ fiR^T, S) 

fl{T,S,t) /i(S') 



(13) 





calculated when all the N{t, S) UPOs are known. 
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Thus, the amount of error that our estimate [Eq. fll3p ] has as when 
compared to the "exact" value of p (when all the UPOs are known) 
can be calculated by 

(16) E[T,S,t\J-^^^ 

which means that the quantity E gives the amount of deviation, in a 
scale from to 1, of p^f [Eq. f[T^ ] as when compared to the "exact" 
value of p [Eq. ((ID])]. Notice that in Eq. (USD, the quantity lOOE 
corresponds to the percentage of error that our estimation has. 

4.3. Uniformly distributed UPOs. There is another way to esti- 
mate the value of p(r, S) in terms of the number of UPOs in a subset 
S" of a chaotic attractor A. We define N{t) as the number of fixed 
points of F'^ in A, N{t, S) as the number of fixed points of F'^ in S, 
Npi{t, S) as the number of fixed points of F"^ in S whose orbit under 
F is recurrent and A^7v/?(t, S) as the number of fixed points of F'^ in 
5* whose orbit under F is non- recurrent. Then, for a sufficiently large 
T and for a uniformly hyperbolic dynamical system for which periodic 
points are uniformly distributed in A, we have 

Using the previous approximations we can write 

_ Nr{t,S) N^r{t,S) ^ N{t,S) 
N{t) N{t) N{t) ■ 

By Eq. ffTOj) we may write p(r, S) ~ Pnr{t, S) and we have that 

(17) p(r,S)«,[S)-^^. 

N[t) 

which can be written as 

Nr{t,S) 
N{r,S) 

Again, we have an expression with a quotient 

Nr{t,S) 
iV(r, S) 

that is, again, easy to obtain numerically by the same technique from 
which pr/ p can be estimated and therefore we can obtain an estimation 
for p(r, S*), represented by pn, by 

(19) p^[r,S,r{t,S)]= p{S)[l-'^^^ 

where r{t, S) represents the number of recurrent UPOs out of a total 
of t UPOs, exactly as previously defined. 



;i8) p{r,S)^p{S)[l 
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5. KOLMOGOROV-SlNAI ENTROPY 

In 1958 Kolmogorov introduced the concept of entropy into ergodic 
theory and this has been the most successful invariant so far[16j. In 
this section we explain how to use the density of first Poincare returns 
to estimate the Kolmogorov- Sinai entropy Hks- 

The exposition here does not aim to be rigorous, only to explain 
how we have arrived at the numerical estimates for the logistic map of 
Sec. El which is a non-uniformly hyperbolic map. 

It is known that|17j 

(20) iV(r) ocexp(ri7^s). 

Consider F as a dynamical system that has the following property: 



N{t) 

for a sufficiently large r. For example, dynamical systems for which 
periodic points are uniformly distributed on the chaotic attractor A 
have this property. 

Considering the tent map and S C [0,1] such that Nnr{t, S) = 1 (if 
there is more that one non-recurrent UPO of period r inside 5" we shrink 
S to have only one), we have p(r, S) ^ ^ that agrees with example l3.H 
for r close to Tmin and for most intervals S. For other non-uniformly 
hyperbolic systems as the logistic the Henon maps, this property holds 
in an approximate sense and this approximation is better the larger r 
is and the closer the interval 5 is to a Markov partition. 

Using the last approximation together with Eq. fl20l) we may write 

p{T, S) 

for some positive constant h & R. So, we have that 
(21) g,„«ilos(tf'^^>)^i|og^^--"'"^>^ '°^"'> 



r V bp{T,S) J T \ p(r, 5) / r 
We define the quantity H{t, S) as 

(22) if(r,5) = llog^^^^^^'^) 



and then, for 6 > 1, it is clear that 

„ 1 , f NKH(T,S) \ 



so H{t, S) is a local upper bound for the approximation of Hks, con- 
sidering a sufficiently large r. 
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Supposing that there is at least one non-recurrent UPO inside S, 

fNRir,S 
Pir,S) 



then for large r we have ^ 6, as 6 is constant. Thus, the term 



1, fNMR{T,S) 

log 



dominates the expression fl2T]) . for longer times. 

This equation allows us to obtain an upper bound for p(r, S). See 
that p(r, S) < Nnr{t, S) exp{-THKs) and if r Tmin then Nnr{t, S) ^ 
1 and we obtain p(r, S) < exp{—THKs) as in example 13.11 

Equation (^2^ depends on the choice of the subset S and is then a 
local estimation for H^s- To have a global estimate we take a finite 
number, n, of subsets Si in the chaotic attractor and make a space 
average as 

'NNii{T,Si[ 



1 

(23) -J]lo, 



Better results are obtained taking the average over pairwise disjoint 
subsets Si that are well distributed over A. 

When we consider A^7v/?(t, S) = 1 this means that we have only one 
non-recurrent UPO, with period r, inside S. In general, for sufficiently 
small subsets. Si, we may have Nnr{t, Si) = 1 Wi and we obtain an 
approximation that only depends on the density function of the first 
Poincare returns 

(24) HKs^—J2^og(-^ 

rn ^ \p{r, S, 

An equation which can be trivially used from the experimental point 
of view since we just need to estimate p(r. Si) and we do not need to 
know the UPOs. For practical purposes, we consider in Eqs. fl22|) . (|23|) 
and that r = Tmin- 

6. Numerical results 

We illustrate our ideas with simulations on the logistic family F : 
[0, 1] ^ [0, 1] given by 

(25) F{x) = cx{l- x), 

were c E R. There are many biological motivations to study this family 
of maps[TS]. The maps that we obtain when the parameter c is varied 
have interesting mathematical properties. It is therefore of relevant use 
for mathematical and biological study. Moreover, for this family it is 
possible to compare the estimates made using all the UPOs to those 
using only some UPOs. 

For most numerical simulations in this section we take c = 4 in Eq. 
(125|) . for which the map is chaotic and the chaotic attractor is compact. 
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Figure 2. Density function of the FPRs, p(r, 5), as 
empty circles and the measure of the non-recurrent pe- 
riodic orbits, /iAr/j(r, 5), as crosses, considering the fol- 
lowing intervals: (A), S = [0.3 - 0.05,0.3 + 0.05]; (B), 
S = [0.3 - 0.01, 0.3 + 0.01]; (C), ^ = [0.3 - 0.005, 0.3 + 
0.005]. 

6.1. Calculating p when all UPOs are known. Figure [2] shows 
the function p(r, S) calculated by Eq. (127|) and the values of UNnij, S) 
calculated by Eq. ([8]), for some subsets S. See that the DFP can be 
almost exactly obtained if all the non-recurrent UPOs inside S with 
period r can be calculated: In Sec. [3] we concluded that p(r, S) ~ 
fJ'Nni'T, S). 

6.2. Calculating p when not all UPOs are known. Figure [3] shows 
the approximations for p{t, S) using Eqs. (|T3l) and (|T9i) . In (B), com- 
paring with (A), we consider longer first return times. We only use 
Eqs. (USD and (J^ for r > 2r„i„. 

6.3. Error of our estimation when not all UPOs are known. 

To numerically calculate the error [Eq. (|T6l) ] of our estimation in Eq. 
(fT3|) . we only consider UPOs with a period smaller than 20. The rea- 
son is because in order to calculate the quotient q2 in Eq. (fT5|) . all 



DT7atc:ttv np ptrqt rtttttrmc: PTTRinnTn orrttc: AMn uo. TTATTRnpv ir 




Figure 3. Red empty circles represent p{t, S) esti- 
mated by Eq. f|T2l) . green crosses estimated by Eq. 
(fTSl) and the black line calculated by Eq. flTTI) . Pic- 
ture (B) is just a similar reproduction of (A) considering 
longer first return times. We consider 200 UPOs inside 
S = [0.1 - 0.001, 0.1 + 0.001], for each r. 



the UPOs must be known. Considering larger periods than 20 would 
be computationally demanding, even thought the proposed method to 
calculate UPOs is capable of finding them all. 

It is also required that r > 2Tmin, once that to calculate the quotient 
qi in Eq. f|T^ there has to exist at least one recurrent UPO within 
the set of t UPOs considered, i.e. r > 1. Therefore, we need to 
choose the size of the interval such that 20-2rmm — 1 is sufficiently 
large, meaning an interval for which Tmin is sufficiently smaller. We 
have chosen e=0.02. 

Since the error of our estimation is proportional to a quotient be- 
tween two quantities that depend on the number r of recurrent UPOs, 
it is advisable that one consider intervals for which a reasonable number 
of recurrent UPOs are found, even when their period is short (smaller 
or equal than 20). Such interval is positioned in places were the nat- 
ural measure is large. In the case of the logistic map, these intervals 
are positioned either close to x=0 or x = 1. Therefore, we consider an 
interval positioned at x = 0.04. From the previous considerations, we 
consider that the interval has a size of e = 0.02. 
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Figure 4. We show the quantity E[t, S, t] with respect 
to the number t of UPOs randomly chosen, for r = 9 
(A), r = 10 (B), r = 11 (C), r = 12 (D), r = 13 (E), 
r = 14 (F), r = 15 (G), r = 16 (H), r = 17 (I), r = 18 
(J), r = 19 (K), and r = 20 (L). The quantity E gives 
the amount of deviation, in a scale from to 1, of pm 
[Eq. (fT3|) ] as compared to the "exact" value of p [Eq. 
(fTOl) ]. We consider an interval positioned at x = 0.04 
with size e = 0.02. 



In Fig. Hl^A-I), we show the quantity E[T,S,t] with respect to the 
number t of UPOs randomly chosen, for r = 9 (A), r = 10 (B), r = 11 
(C), r = 12 (D), r = 13 (E), r = 14 (F), r = 15 (G), r = 16 (H), 
r = 17 (I), r = 18 (J), r = 19 (K), and r = 20 (L). 

The most important information from these figures is that as UPOs 
of longer periods are considered [going from Fig. (A) to (L)], the error 
E of our estimation decreases in an average sense considering all the 
values of t. Another relevant point is that the larger the number t of 
UPOs considered, the smaller the error. Notice that the total number 
of UPOs of period r is given by 2^. Therefore, looking at Fig. lU^L), 
one can see that even considering only about 0.0009% of all the UPOs 
(10 UPOs, out of a total of 2^^=1048576), the error of our estimation 
is smaller than 14% when compared to the "exact" value of p. 




Figure 5. (A) A bifurcation diagram as points (light 
green) and the randomly chosen intervals as empty 
(black) squares. (B) Lyapunov exponent as line and filled 
circles representing the H^s entropy using Eq. fl22]) . for 
the logistic family. We consider 400 values of c and for 
each c the size of the set S" is e = 0.002. 



6.4. Estimating the KS entropy. In order to know how good our 
estimation for Hks is we use Pesin's equahty which states that Hks 
equals the sum of the positive Lyapunov exponents, here denoted by A. 
For the logistic map there is at most one positive Lyapunov exponent. 

Figure [5] shows the approximation for the quantity Hks using Eq. 
fl22|) . See that Eq. fl22l) only needs one subset S on the chaotic attractor 
to produce reasonable results. In this numerical simulation we vary the 
parameter c of the logistic family and for each c we use just one subset 
S{c) randomly chosen [shown in Fig. O (A)] but satisfying Tmin G 
[10, 14] so that T considered in Eq. (l22i) is sufficiently large. 

Finally, Fig. [6] shows the global estimation for Hks, using the Eqs. 
(123|) and (124|) . considering 40 intervals Si for each value of c. Recall 
that if A < 0, then Hks = 0. 

6.5. Numerical work to find UPOs. The analytical calculation of 
periodic orbits of a map is a difficult task. Even for the logistic map it 
is very difficult to calculate periodic orbits with a period as low as as 
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Figure 6. The Lyapunov exponent A as line and the 
aproximation of Hks entropy using Eqs. and 
as empty circles. (A), Eq. (B), Eq. (jM]). In this 

simulation we consider 100 values of c and for each c we 
consider 40 subsets Si each one with lenght e = 0.002. A 
subset Si is picked only if r^m £ [10, 14]. 



four or five. In our numerical work we need to find unstable periodic 
orbits and, in some cases, we need to find all different UPOs inside a 
subset of the phase space, for a sufficiently large period. For that, we 
use the method developed by Biham and Wenzel[19j. They suggest a 
way to obtain UPOs of a dynamical system with dimension D using 
a Hamiltonian, associated to the map, with dimension ND, where 
is the number of UPOs with period p. The extremal configurations of 
this Hamiltonian are the UPOs of the map. The force dH/dt directs 
trajectories of the Hamiltonian to the position of a UPO. 

The Hamiltonian associated with the map gives a physical interpre- 
tation of the problem but in some cases it is impossible to know it. We 
propose a method with a similar interpretation that is simpler in the 
sense that we do not need to know the Hamiltonian associated with the 
map, just an array of A^ coupled systems where the linear coupling be- 
tween nodes acts as the force directing the network to possible periodic 
solutions of the dynamical system concerned. 
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For this method we just need the force associated with the ith node, 
described by x\ and satisfying the Euler-Lagrange (E-L) equations: 

d dL _ dL 
dt dx^ dx^ ' 

where L is the Lagrangian associated with the map. We are inter- 
ested only in static extremum configurations of the Hamiltonian and 
therefore the kinetic term will be neglected This implies 

— - 

dx^ 

We illustrate the numerical calculation of UPOs with arbitrary length 
applying it to the logistic family. Because the static (E-L) equations 
reproduce the map, we have 

dL 

'Q~i ~ -^n ~ — x^. 

The force of the i node will be given by 

dL ■ ^ 

When the chain is in stable or unstable equilibrium (an extremum 
static configuration of the Hamiltonian), Fi = for all i. To find 
a specific extremum configuration of order p of the Hamiltonian we 
introduce an artificial dynamical system defined by 

dx^ 

(26) ~dt^''-^'' ^ = 1,-,P, 

where Sj = ±1 represents the direction of the force with respect to 
the ith node. This equation is solved subject to the periodic bound- 
ary condition x"^^^ = x^ and when the forces in all nodes decrease to 
zero the resulting structure x^ is simultaneously an extremum static 
configuration and an exact p-periodic orbit of the logistic map. For 
c = 4, if we take Sj = — 1 Vz then we obtain the trivial periodic point 
Xi = Vi. The different ways to write Si will give different UPOs. We 
may look at Sj as the representation of the orbit in a symbolic dynamics 
with S = { — 1, 1}, taking the trivial partition on the logistic map, i.e.. 
Si = -1 if Xi e [0, 1/2] and = 1 if G [1/2, 1]. 

Equation fl2B]) is in fact an equation for a network of coupled maps. 
The UPOs with period p embedded in the chaotic attractor can be 
calculated by finding the stable periodic orbits of the following array 
of maps constructed with i = l,...,p nodes x^, where every node is 
connected to its nearest neighbor as in 

■^n+l ~ -^n ~ CSj[x„ — -F(x„)], 

with the periodic boundary condition x^ = x^, where the term csi[xl^^ — 
F{xl^)] represents the Lagrangian force. 
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7. Conclusions 

In this work we propose two ways to compute the density func- 
tion of the first Poincare returns (DFP), using unstable periodic orbits 
(UPOs), where the first Poincare return (FPR) is the sequence of time 
intervals that a trajectory takes to make two consecutive returns to a 
specific region. In the first way, the DFP can be exactly calculated 
considering all UPOs of a given low period. In the second way, the 
DFP is estimated considering only a few UPOs. We have numerically 
shown that the error of our estimation becomes smaller, the longer the 
period of the UPOs and the larger the number of UPOs considered. 

The relation between DFP and UPOs allows us to compute easily 
an important invariant quantity, the Kolmogorov- Sinai entropy. 

For non-uniformly hyperbolic systems there exists some particular 
subsets for which the UPOs that visit it are not sufficient to calculate 
their measure of the chaotic attractor inside it [201 [21] . For such cases 
our approach works in an approximate sense, but it still provides good 
estimates as we have shown in our simulations performed in the logistic 
map, a non-uniformly hyperbolic system. In addition, the approaches 
shown in here were applied in ref. |22] to estimate the value of the Lya- 
punov exponent in the experimental Chua's circuit and in the Henon 
map, both systems being non-hyperbolic. 

Our approach offers an easy way to estimate the KS entropy in ex- 
periments, since one does not need to calculate UPOs, but rather only 
to measure the DFP of trajectories that make shortest returns, i.e. the 
quantity p{Tmin, S). These are the most frequent trajectories, and as a 
consequence even if only a few returns are measured, one can obtain a 
good estimation of p^Tmm, S). More details of how to estimate the KS 
entropy from experimental data can be found in Ref. f27\ . 

8. Appendix 

8.1. Measure and density in terms of FPRs. We calculate p(r, S) 
also in terms of a finite set of FPRs by 

(27) ^l--^)-^ 

where K{t, S) is the number of FPRs with a particular length r that 
occurred in region 5* and L{S) is the total number of FPRs measured 
in S with any possible length. 

We calculate fi{S) also in terms of FPRs by 

(28) piS) = ^ 

riL 

where ul is the number of iterations considered to measure the L{S) 
FPRs and so til = J2t=i '^n (see definition 13. ip . 
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We define the average of tlie returns by 

Comparing Eqs. (I25|) and f l2I?]) . we liave tliat 

(30) ^^{S) = -i- 

< r > 

also known as Kac's lemma. 
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